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ABSTRACT 

Metallic gas detected in the Pictoris circumstellar debris disk raises many questions. The origin of this 
gas is unclear and its very presence is difficult to explain: many constituents of the gas are expected to be 
radiatively accelerated outward, yet their motion appears to be consistent with Keplerian rotation out to at least 
300 AU. Hydrogen has previously been hypothesized to exist in the disk, acting as a braking agent, but the 
amount required to brake individual elements conflicted with observed upper limits. 

To resolve this discrepancy, we search for alternative braking mechanisms for the metallic gas. We find that 
all species affected by radiation force are heavily ionized. Frequent Coulomb collisions couple the ions into 
a single fluid, reducing the radiation force on species feeling the strongest acceleration. For a gas of solar 
composition, the resulting total radiation force still exceeds gravity, while a gas of enhanced carbon abundance 
could be self-braking. We also explore two other braking agents: collisions with dust grains and neutral gas. 
Grains surrounding (3 Pic are photoelectric ally charged to a positive electrostatic potential. If a significant 
fraction of the dust grains are carbonaceous (10% in the midplane and larger at higher altitudes), ions can be 
slowed down to satisfy the observed velocity constraints. In this case, both the gas kinematics and spatial 
distributions are expected to coincide with those of small grains, the latter being indeed observed. For neutral 
gas to brake the coupled ion fluid, we find the minimum required mass to be w 0.03 M®, consistent with 
observed upper limits of the hydrogen column density, and substantially reduced relative to previous estimates. 

Our results favor a scenario in which metallic gas is generated by grain evaporation in the disk, perhaps 
during grain-grain collisions. We exclude a primordial origin for the gas, but cannot rule out the possibility 
of its production by falling evaporating bodies near the star. We discuss the implications of this work for 
observations of gas in other debris disks. 

Subject headings: stars: individual ((3 Pictoris) — circumstellar matter — acceleration of particles — scattering 
— planetary systems: formation — planetary systems: protoplanetary disks 



1. INTRODUCTION 

The evolution of circumstellar disks is intima tely connected to the pro cess of star and planet formation. In the standard 
isolated low-mass star formation scenario (e.g., Shu, Adams & Lizano 1987), disks begin their existence as protostellar accretion 
disks, later turning into protoplanetary disks, once the newborn star reaches the T Tauri phase. The late evolutionary stages 
of circumstellar disks are not well understood. Once the star has reached the main sequence and the planet formation process 
is complete, the disk is thought to have been depl eted of its prim ordial gas, keeping small amounts generated by either stellar 
winds or evaporation of solid bodies (e.g., [Lagrange et al. 2000). In the few cases where late disks are observable, lifetimes 
of small dust grains (< 1 /i m) are shorter than the estimated ages of those systems, implying a dust replenishment mechanism 
(Backman & Paresce 1993). These kind of late circumstellar disks are usually referred to as debris disks. The gas lifetime in 
those systems, currently unknown, determines the maximum period over which protoplanetary cores can accrete gas and form 
gaseous planets (e.g., Bodenheimer & Lin 2002), and is related the formation of terrestrial planets and cores of giant planets by 
the increased efficiency of planetesimal accretion due to the presence of gas (lRafikovl20 04| ). 

So far, the best studied debris disk is that of (3 Pictoris, a nearby (19.3 ± 0.2 pc. ICrifo et al]ll997l) . A5V star. The dust 
comp onent of the disk has been extensively studied since it revealed itself as infrared excess in the stella r spectrum lAumann et alJ 
1985). The fav orable edge-on orientation has allowed d irect imaging in thermal emission (Lagage & Pantir3 ll994l) . as well as in 
scattered light ( Smi th & Terrilelll984l iHeap et alJl2000l) . indicating a disk size ~1 000 AU. The gas component has been more 
elusive. Freudling et al. ( 1995) set a 3 a u pper limi t on the column d ensity of H I, A/h i < 10 19 gem" 2 , by non-detection of the 21 cm 
line, while Thi et al. (2001) and Lecavelier des Etangs et al. (2001) searched for signatures of H2 in the IR and UV, respectively, 
obtaining discrepant results: 50M ffi in the first case (corresponding to Nh, ~3x 10 21 cm" 2 ), and an upper limit Nh^ < 10 18 cm~ 2 
(~0. 1 Mm) in the second. Metallic elements have been detected through stable and variable circumstellar lines ( Hobbs et al.ll985t 
iFerlet et al]|1987l) . The variable, redshifted lines have been attributed to falling evaporating bodies (FEBs) in high ly eccentric 
orbits (Ferlet et al. 1987), whereas the stable component has no well established explanation for its origin. Lagrange etaLl(Ti998l) 
proposed three possible gas formation scenarios: stellar wind, star-grazing comet evaporation, and dust grain evaporation near 
the star. Recent improvements in the gas characterization inc lude spatially resolve d emission lines from stable metallic species, 
which revealed the m to be con sistent with Keplerian rotation (Olofss on et al.l20 01). and spatially extended out to at least 300 AU 
from the star ( Brande ker et all2004|) . 

It was realized early on (e.g.. lBeust et al.ll9 89) that certain elements in the gas disk of (3 Pic should be subject to a very strong 



e-mail: fernandez@ast.ro. utoronto.ca 



2 



Fernandez, Brandeker, & Wu 



radiation force from the star, overcoming gravitational attractio n. However, b esides Keplerian rotation, stable lines do not show 
any substantial radial velocities relative to the star (e.g., Lagran ge et all 19 98. Brandeke r et all2004l) . This implies the presence 
of a braking mechanism acting in the disk. 

The first model of how gas might be braked and produce the stable absorption lines was presented by Lagra nge et alJ (119981) . 
They proposed a scenario in which gas is injected into the disk due to FEBs. Since neutral hydrogen is not affected by the stellar 
radiation force, it would accumulate in the region of gas generation, in the form of an annulus located a few AU from the star. The 
elements affected by radiation force would then be slowed down temporarily in this annular region, producing the observed stable 
absorption features. After passing this hydrogen-rich region, elements would then be accelerated again, leaving the system at 
high velocities. This model reproduces the stable absorption lines, and is consistent with the FEBs scenario for the variable lines. 
However, the fact that elements affected by radiation force reach high velocities after passage by the ann ulus, which extends only 
out to a few AU from the star, is incompatible with the emission from ~300 AU. Brandeker et al. (2004) calculated the amount of 
material necessary to slow down the most strongly accelerated particles to with in observational constraints b y means of neutral- 
neutral collisions, finding that ~50M ffi of hydrogen are needed. According to Thebault & Augereau (2005), there could be no 
more than 0.4M ffi of gas in the disk, as it would affect the dynamics of the dust to the degree that dust production models become 
incompatible with observations. 

As an alternative braking mechanism. lBrandekerl (120041) explored the slowing down of ionized s pecies by a large-sc ale magnetic 
field, concluding that only a toroidal field could remove the radial velocity components of ions. [Brandeker (2004) showed that 
unless the field is stronger than ~1 mG at 100 AU, this configuration would be unstable to currents generated by radiation force. 
He also raised the question of the origin of such an ad-hoc field geometry. 

In this paper we address the braking of the stable metallic gas in the entire Pic disk. Given the spectral type of the star and the 
fact that the disk is optically thin, we expect several gaseous species to be ionized. This prompts us to explore the role of Coulomb 
collisions among ionized gas particles, which may be an important braking agent due to the long range of electromagnetic 
interactions. We proceed then to explore the interaction of this ionized gas with dust grains, which are photoelectrically charged. 
Finally, we extend the study of Lagrang e et al.l 1199 8) to explore whether ion-neutral collisions are able to brake particles affected 
by radiation force far out in the disk, and the amount of neutral gas required to satisfy the observed velocity constraints. 

The paper is organized as follows: in ^2] we explore the physical conditions of the gas, and the constraints on the possible 
braking mechanisms that this imposes. In 5j5]we study ion-ion collisions, ion-grain collisions, and revisit ion-neutral interactions. 
In ^|we discuss the implications of our results for gas generation and extend it to similar systems. Our conclusions follow in |5] 

2. SETTING THE STAGE FOR BRAKING 

In this section we explore the physical environment of the gaseous disk. This will enable us to characterize the requirements 
the braking mechanisms have to satisfy. 



where c is the speed of light, is the cross section for the transition between atomic levels j and k integrated over wavelength, 
and F\ is the stellar flux per unit wavelength evaluated at the line center. Since the stellar flux decays like the inverse square of 
the distance from the source, the ratio (3 between radiation and gravitational forces is distance-independent: 



where we have expressed ajt in terms of known atomic quantities (e.g.. lHilbornll982T) . Here, gj and gu are the statistical weights 
of levels j and k, A^j is the Einstein coefficient for spontaneous emission from level k to j, \jk the wavelength of the corresponding 
transition, G is the gravitational constant, M the stellar mass, and m the mass of the atomic species under consideration. 

The energy of a particle of mass m in a bound orbit with semi-major axis a around an object of mass M is £bind = —GMm/(2a). 
When the contribution of radiation is added, the energy of the system increases by an amount E m & = j3GMm/r, where r is the 
distance to the star at a given point in the orbit. For a circular orbit (r = a), the particle ceases to be bound if (3 > 0.5. If (3 < 0.5 
the particle remains in a bound orbit, but "feeling" a lower dynamical mass (1 - /3)M. 

Values for j3 in this system are available in the literature for some elements, but we find it necessary to calculate a set of (3 
values for as man y elements as possible. The low temperatures and densities of the gas in the disk (at distances > few AU from 
the star, see ^2.21 imply that radiative de-excitations are faster than either collisions or radiative excitations, thus most atoms are 
in the ground state. This impli es j = in equation (0, simplifying calculations considerably. Since we do not expect species to be 
more than twice ionized (see 32. 2> . we study the first three ionization states of elements from H to Ni. The atomic data are taken 
from the NIST Atomic Spectra Database version 2.0 (Mar tin et al.ll999l) 2 . For element s with multiple groun d states (e.g., Fe I), 
we assume that the population of each level is proportional to its statistical weight (e.g., Lagrange et al Jl996l) . The uncertainties 
in A*o range from <3% for the strongest lines to ^100% for the weakest 

The stellar flux is approximated with a PHOENIX model spectrum ( H auschildt. Allard & BarorJll999 ) for a star of eff ective 
temperature T e g = 8000K and surface gravity logg = 4.2. Given that /3Pic is rotating at 130 ± 4km s -1 ( Rover et afll2002l) . the 



2.1. Radiation Force 



The radiation force F m & acting on a given atom originates by a net momentum transfer due to scattering events, 




(1) 




(2) 



2 http://physics.nist.gov/asd2 
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TABLE 1 

Radiation force coefficients obtained from rotationally broadened /3 Pic-like spectrum 



Ion 


(3*- b 


Ion 


/3 b 


Ion 




Ion 


F» 


HI 


(1.6±0.1)10 -3 


FI 





SI 


0.56 ±0.09 


VI 


72±4 


He I 





FII 


(3.5±0.9)10~ 6 


SII 


(9.0±1.0)10~ 5 


VII 


4.4 ±0.2 


He II 




Fill 


(5.0±1.0)10 -9 


S III 


(2.0 ±1. 0)10^ 


VIII 





Li I 


900 ±40 


Nel 





CII 


(2.3±0.4)10" 3 


CrI 


93 ±5 


Li II 





Ne II 





CI II 


(3.7±0.4)10~ 7 


Crll 


(6.0±3.0)10~ 7 


Li III 




Ne III 


(9.0±2.0)10~ 8 


CI III 


(3.0±2.0)10~ 6 


Ci'III 




Bel 


62 ±7 


Nal 


360 ±20 


Arl 


(1.7±0.3)10~ 6 


Mill 


28 ±3 


Bell 


124 ±6 


Na II 





Aril 





Mn II 


7 ± 1 


Belli 





Na III 





Arlll 


(1.5 ±0.2)10~ 7 


Mnlll 




tJ T 
D 1 


on 1 1ft 
JU ± IU 


Mgl 


HA _1_ Q 


1/ T 
IS. 1 


TAA _L TA 


re 1 




BII 


0.07 ±0.04 


Mg II 


9±2 


KII 




Fell 


5.0 ±0.3 


B III 


19± 1 


Mg III 


o 


Kill 


(A a + n 21 1 0~ 4 


Felll 


f-j n + o 61 1 0~ 7 


CI 


(3.3±0.1)10 -2 


All 


53±6 


Cal 


330 ±40 


Col 


16±1 


CII 


(2.3±0.2)10~ 3 


Al II 


0.36 ±0.05 


Call 


50±10 


Coll 





cm 


(8.5±0.9)10~* 


Aim 


12±1 


Calll 




Colli 


(4.0±2)10- 7 


NI 


(2.1±0.1)10 -4 


Si I 


6.0 ±0.6 


Scl 


220 ± 20 


Nil 


26±2 


Nil 


(7.5 ±0.5)10"* 


Sill 


9±9 


ScH 


(1.3±0.4)10 3 


Nill 


(7.0±2.0)10~ 2 


NIII 


(7.0 ±1.0)10-* 


Si III 


(5.8±0.6)10~ 4 


ScIII 


(9.0±3.0)10~ 2 


Nilll 


(3.0±2.0)10~ 7 


OI 


(3.3±0.2)10" 4 


PI 


3.4 ±0.6 


Til 


97±5 


P° 


4.4 x 10-" 


on 


(3.1±0.7)10" 9 


PII 


(2.2 ±0.3) 10~ 3 


Till 


28±2 


e c 


0.27 


OIII 


(6.5±0.6)10~ 7 


PHI 


(5.0±2.0)10~ 4 


Tilll 


(5.0±0.1)10~ 4 







11 /3 = means that no ground state transitions are known in the range 1 000 A < A < 50 000 A 
b Empty entries mean that no atomic data are available 

c Proton (H II) and electron values are calculated using Thomson cross section 



spectral features are significantly broadened. A proper rotation of the spectrum requires detailed resolution of the photospheric 
absorption profiles, thus a special version of the model atmosphere was used in this work, which has constant spacing in wave- 
length, namely 0.05 A in the region 1 000-10 000 A, and 0.5 A in the region 10 000-5 000 A (P. H auschildt & I. Kamp, private 
communication). A rotational broadening of this model spectrum was performed (e.g.. lGravfl 976). assuming a linear limb dark- 
ening law with e = 0.5. This choice adds an uncertainty in the relative flux level at the bottom of photospheric lines ranging 
from 2% in the visual to 10% in the UV. The spectrum was then flux calibrated using photometric information from the Tycho-2 
Catalog (Hdg et al. 2000) and the distance to the star. We estimate the error of this flux calibration to be ~4%. 

Tab leQJshows the results of our calculations for elements from H to Ni. The mass of the star used is M = 1 .75 M JCrifo et alJ 
QUI. Calculation of uncertain ties is exp l ained in Appendix|A| The agreement level with previously calculated values is highly 
variable. As an example, Lagrange et al. ( 1996) calculated /3 values using calibrated HST spectra. Our results agree with their 
values for Fell, A1II and Al III within our uncertainties, but their result for Mgll differs from ours by a factor of 2. We take 
the good agreement between their Fe II value and ours as a consistency check, given the great number of transitions required 
to calculate (3 for this ion. Nevertheless, when looking at our results, it should be kept in mind that they were obtained with a 
generic model spectrum, therefore subject to possible systematic deviations from the particular [3 Pic case. 

2.2. Ionization State of the Gas 

Since we are interested in relating the dynamics of the gas to electrostatic interactions, we need to analyze in detail the 
ionization state of the gas. Given that the disk is optically thin, and the densities are low, the ionizat ion state is de termined by 
photoionization and radiative recombination. The temperature of the gas is below 300 K (Kamp & van Zadelhoff 2001), thus 
recombination rates are not very sensitive to the exact value of the temperature, depending most strongly on the electron density. 
Given the spectral type of the star (A5V), we expect most elements to be either singly ionized, or neutral. 

We thus wrote a simple photoionization code which calculates the densities of neutral and ionized elements from H to Ni 
for an optically thin gas disk, with the corresponding electron density, and assuming solar composi tion. The latter cho ice is 
motivated by observations showing refractory elements having solar abundance relative to each other ( Laaranae et al. 1995J). The 
stellar ionizing radiation is obtained from our rotated and flux-calibrated mod el spectrum, whereas i nterstellar bac kground UV 
field and and cosmic ray ionization rate (2 x 10~ 17 s _1 atom" 1 ) were taken from Weingartner & Draine (2001) and Spitzer ( 1978), 
respectively. Radiative recombination coefficients and photoionization cross sections were taken from the Cloudy C96.01 atomic 
data library ( Ferland et al. 1998). Recombination of ions on dust grain surfaces, which can be the dominant recombination process 
in the interstellar medium, is not included as the dust in the disk is expected to have a highly positive charge (Appendix 151. in 
contrast to the negatively charged grains of the interstellar medium. 

At low temperatures the chemistry of the gas is not affected significantly if one sets T %3S = 7d ust ( iKamp & van Zadelhoffl200"l . 
where r gas is the gas temperature and rd ust is the dust temperature. We thus adopt a temperature profile corresponding to grains 
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p[AU] 



FIG. 1. — Plotted against cylindrical radius p are: average midplane dust temperature T (eq.|4] right axis), which we take t o be the same as the local 
gas temperature; and midplane number densities (left axis) of electrons and selected species. Nal is obtained from observations i Brandeker et al. 2004), the 
remaining species are obtained from our photoionization code. C is partially ionized (~50%) throughout the disk (and the main contributor of electrons), while 
Na and Fe are almost completely ionized. 



of size a (e.g., Kamp & van Zadel hoffl2001l) 



r a = 282.5 ( — 



1/5 



AU 
r 



2/5 i 



(3) 



where L is the stellar luminosity (11 L Q for /3Pic), and r is the distance to the star. To eliminate the dependence on grain size, we 
take a mean temperature weighted by grain surface area and size distribution 



J T^a 2 dn 
f a 2 dria 



(4) 



We adopt dn a cx a -3 5 and perform the integral over the range a m ; n < a < a raax , with a^n = 1 Atm and a max = 1 km. The chosen size 
distribution implies that the total grain surface area, and therefore our weighted temperature, is dominated by values correspond- 
ing to the smallest grains, being insensitive to the choice of a max . The resulting temperature profile in the disk midplane is shown 
in Figure [2 (right axis). 

The observational constraint on the gas density comes from the fit to Nal emission (Br andeker et alJ2 004), which has the form 



«NaiO,z)= 10 



0.94 



6.32' 



"1-1/2 




exp 


_ (o.l7p) _ 



(5) 



where p is the distance to the star in the disk midplane, z the height above the midplane (r 2 = p 2 +z 2 ), and po = 117 AU. The Na I 
profile in the disk midplane is shown in Figure ^ scaled by a factor 10 5 . Our photoionization code numerically solved for the 
densities of neutral and ionized species that are consistent with this Na I profile, the temperature profile of equati on l|4 ), and our 
assumption of solar composition. The sensitivity of our results to the adopted temperature profile is discussed in 34.11 

Our results show that H, He, N, O, F, Ne, and Ar, if present, are almost completely neutral, their marginal ionization (^10~ 6 ) 
being determined by interstellar UV flux and cosmic rays. This group is unaffected by radiation force in its neutral phase. The 
elements C, CI, Be, P, and S are partially ionized, with ionization fractions ranging from ~50% for carbon to >95% for beryllium. 
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TABLE 2 

Velocities before ionization for neutral 
species with p > 0.5 



Element r sta r a 
(s- 1 ) 



Li I 


9.0 x 10~ 


-6 


Be I 


1 1 x 10 


-9 


B I 


7 5 x 10" 


-8 


Nal d 


1.1 x 10" 


-7 


IVTit T 
Mgl 


f\ 8 v 1 fl" 


-8 


All 


1.1 x 10" 


-4 


Si I 


3.9 x 10" 


-7 


PI 


4.0 x 10" 


10 


SI 


4.0 x 10~ 


10 


KI 


4.4 x 10 


-7 


Cal 


1.3 x 10- 


-5 


Scl 


6.4 x 10- 


-8 


Til 


1.5 x 10- 


-7 


VI 


4.0 x 10- 


-7 


CrI 


1.0 x 10- 


-7 


Mnl 


6.8 x 10- 


-8 


Fel d 


5.8 x 10- 


-8 


Col 


2.0 x 10- 


-8 


Nil d 


1.0 x 10- 


-7 



Vion b 


ni/«n+ c 


(kms -1 ) 




0.11 


4.4 x 10" 6 


58 


4.4 x 10~ 2 


0.41 


6.1 x 10~ 4 


3.3 


3.2 x nr 4 


1.1 


6.4 x 10~ 4 


5.0 x 10^ 


8.5 x 10" 7 


0.02 


1.3 x 10~ 4 


8.7 


1.5 x 10~' 


1.5 


6.4 x 10" 2 


0.5 


1.5 x 10~ 4 


0.03 


3.5 x 10~ 6 


3.5 


7.0 x 10~ 4 


0.6 


3.0 x 10~ 4 


0.2 


1.2 x 10~ 4 


0.9 


4.8 x ur 4 


0.4 


7.6 x 10~ 4 


0.5 


9.0 x 10~ 4 


0.8 


2.9 x 10" 3 


0.3 


4.8 x 10 -4 



"Stellar ionization rate at 100 AU from the star. 
b Terminal velocity before ionization, independent of 
distance. 

c Ratio of neutral (I) to ionized phase (II or higher) at 
r = 100 AU in the disk midplane. 

d Species with mea sured radia l velocity by 
Brandeker et al. 1 2004 ) [see H2.2I for values and 
uncertainties] 



Inside ~50 AU the stellar ionization dominates, whereas UV background takes over at larger distances. While neutral carbon, 
chlorine, and sulfur feel little or no radiation pressure, the neutral phases of phosphorus and beryllium are accelerated. The 
potential of these species to trace the velocity of neutrals will be discussed below. The remaining species (including Na) are 
ionized to fractions >99.9%. Electron densities in the disk midplane are shown in Figure^ together with density profiles of C II 
and Fe II, which are the most abundant species by mass in the partially-, and heavily-ionized group. 

The group with ionization >99.9% contains all species affected strongly by radiation force. This is crucial for explaining 
the dynamics of the gas. From Tabled we see that for most of these species, the highest (3 corresponds to the neutral phase 
(exceptions being Be and Sc), which has a very short lifetime compared to that of the ionized phase. Once ionized, elements 
feel a lower radiation fo rce a nd can undergo Coulomb interactions with other charged particles and lose the excess momentum 
acquired from photons ( SI3.1I >. 

A quantity of relevance is the velocity boost a neutral particle acquires before being ionized, v; on . In the optically thin case, the 
ratio between the number per unit time of scattering photons to that of ionizing stellar photons depends only on the shape of the 
stellar spectrum, being independent of the distance to the star. The number of ionizing stellar photons hitting an atom per unit 
time is the stellar ionization rate 

f°° F\ 

r*= / - 7 a ion d\, (6) 
Jo (hc/X) 

where <7i on is the ionization cross section, the ionization time (or lifetime as neutral) being the inverse of r». The velocity boost 
before ionization is then 

GM 1 

Vion^P^- — , (7) 
T 1 * 

where T* is evaluated at distance r from the star. Table|2]shows the results for elements with (3 > 0.5, for which the ionization 
rate due to interstellar UV is negligible relative to r». The most relevant feature is that, with the exception of Bel and PI, all 
neutrals feeling radiation force achieve velocities below 3.5 km s" 1 before being ionized. Since the acceleration in the neutral 
phase is nearly constant, the mean radial velocity is v; on /2, provided no braking mechanism acts on particles in the neutral 
phase. T his is consistent with o bserved values for neutrals: VN a i = 1.2±0.3kms _1 , vpei = 0±0.3kms -1 , and Vnu = 0.4 ± 
0.4 km s" 1 ( Brandeke r et alJl2004h . Hence, as long as ionized particles are effectively braked, there is no need to find a braking 
mechanism for high-/3 neutrals. Observed radial velocities of Be I and P I would be able to test these predictions. Unfortunately, 
their strongest ground state transitions are in the UV, which requires observations from space. The expected line equivalent 
widths, derived using our ionization model, would be 10 mA for P I A 1774, while Be I A2350 would be hardly detectable at 3 jiA. 
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3. BRAKING THE IONIC FLUID 

We have established in the last section that the braking of neutral species is not required, only ions need to be slowed down. 
This changes the nature of the puzzle and makes it solvable. In this section we explore collisional processes involving electrostatic 
interactions that can account for the braking of ionized particles. From simple dynamical considerations it is clear that any braking 
medium needs to satisfy two requirements. First, a high collision frequency with ions, so that momentum can be effectively 
exchanged with particles accelerated by radiation. Second, the medium must provide a high inertia, so that it is not dragged 
along. We thus focus on three processes likely to operate in the disk: ion-ion, ion-dust, and ion-neutral collisions. 



3.1. Collisions between Ions 

The fact that all the species affected by radiation pressure are highly ionized prompted us to explore first the role of Coulomb 
collisions among them. For treating the problem, we assume that magnetic fields are weak enough in the regions of interest, 
~100AU from the star, so that the dynamics of charged particles is not influenced by them. Otherwise , complicated MHD 
effects may come into play. The sta tistical treatment of collisions among charged particles was studied by Spitzer ( 1956), based 
on results of Chandrasek harl (I 1 94 111 1 9431) for gravitational interactions. In Spitzer's framework, particles are categorized as test 
and field particles, the former having an initial velocity relative to the center of mass of the latter. Field particles are assumed 
to follow a Maxwellian velocity distribution. Braking of test particles is then characterized by the time it takes for the beam of 
test par t icles t o diffuse in velocity space. This time corresponds to dynamical friction in the gravitational case, and is given by 
ISpitzerU 19561) as 

t u m?(k B T f )v 

S ^ ' ' 8™ f Z f 2 Z t V( OTt + OTf )//(£)i n A' ^ ; 

where the subscripts t and f correspond to test and field particles, respectively. Here, v is the velocity of test particles relative 
to the center of mass of field particles, m is the mass of individual particles, n is the number density, Z is the charge in units of 
the electron charge e (throughout this work we assume Z t = Zf = 1 for ions), 7f is the temperature of the field particles, and ks is 
Boltzmann's constant. The quantities £ and A are defined by 

v 

£ = (9) 

Cth.f 



2kBT{ rim 
c±,f=\ (10) 



Ad 3 f k\T^ 



Omin 2ZfZ t e j V 7T« e J 

where c t h.f is the thermal velocity of field particles, Ad is the electron Debye screening length, f> m ; n is the minimum impact pa- 
rameter for Coulomb collisions, and T e is the electron temperature. In this context, electrons and ions are sufficiently thermalized 
to allow setting T e = Tf. The In A factor accounts for the many weak scatterings out to distance Ad- The functions H and $ are 
given by 3 

V*" Jo 

Two related questions are considered here. First, for a newly ionized particle with speed v; on , how long does it take to equalize 
its speed with the rest of the ions? Second, are high-/? ions dynamically coupled to low-/? ions? 

The first question can be answered by evaluating f§ with the relative velocity v = vi on , where v; on is the velocity of neutral 
particles just before ionization, shown in Table|2] The resulting fs in the disk midplane is shown in Figure|2]in units of the local 
Keplerian orbital period. It takes a very short amount of time for most newly ionized particles to be equalized with the ionic sea 
through Coulomb collisions, except for the element Be. A prediction is therefore that Be will be preferentially lost from the disk. 

To address the second question, we express the timescale to radiatively accelerate a test ion as 

_ «tV« _ r 2 v* 

acc F md (A-l)GM' (i4) 

where v* is the radial velocity relative to the star, and f3 t > 1 is the test particle radiation force coefficient. If fs *C face, any excess 
momentum acquired by species that feel the radiation force will be transferred to the field particles as a result of many random 
collisions, the ions thus being dynamically coupled despite their different (3 values. Notice that the ratio fs/f a cc does not depend 
on v except through the function //(£ = v/cth,f), which reaches a maximum when £ ~ 1. In Figure|5] we plot the ratio fs/?D as a 
function of the disk cylindrical coordinates p and z, taking f = 1, with C II as field particles and Fe II as test particles 4 . The species 

3 In i Spitzer 1956) the function in equation 1121 is called G, but we rename it H to avoid confusion with the gravitational constant. 

4 When ( < 1 or ^ > 1, <s can become longer that r acc . In the former case, ions may drift relative to each other, but with a speed orders of magnitude below 
the local sound speed. 
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FIG. 2. — Ratio of slow-down to orbital timescales along the disk midplane, for ions with incident velocities v = Vj n, where v; on is the velocity of neutral 
particles just before ionization, shown in Tablel2l We adopt C II as field particles, but the results remain similar if other field particles are used. 



chosen are representative as they are the most abundant (in number and mass) within groups that feel or do not feel radiation 

force (/?cn ~ and /3p e n ~ 5, respectively). 

Of particular interest is the coupling of ions at high altitudes over the midplane. Brandeker et al. ( 2004) detected Ca II emission 
being consistent with Keplerian orbits even at a height of z ^80 AU (p ^120 AU). In that zone we have ?s/f a cc ~ 0.05 for v = Cth,f, 
i.e., ions are still coupled there. The effect of radiation force on the coupled ensemble of ions can be quantified by an effective 
radiation force coefficient 



&ff = 



gravj 



(15) 



where F rac y and f gm v.j are the radiative and gravitational forces acting on species /' within a fluid element, (3j is the radiation 
pressure coefficient for species j, and gj is the mass density of species j. Thus /3 e ff is a weighted average of the values for each 
particle, the weight being the mass density of each species. If /3 e ff < 0.5, ions brake by themselves. Since ions outweight electrons 
by a factor >2000, and there is one electron per ion, we neglect the electrons in the estimate of f3 e g and assume that they are 
dragged along the ion fluid, preserving charge neutrality. 

For solar composition, the most abundant ions by mass are, in decreasing order, C(/3 ~ 0), Fe(5), Si(9), Mg(9), S(0), Ni(0) 
and Ca(50). These particles dominate the value of /j e ff, with carbon and iron being the relevant species. If all of them were 
100% ionized, we would have (3 e s w 3.5, but since the ionization fraction of CII varies with distance from the star, the actual 
value of /3 e ff is higher (by a factor 2 at most) and has a weak dependence on position. In Figure|4]we plot /3 e ff averaged over all 
ionized species as a function of disk cylindrical coordinates p and z- A straightforward way of braking the gas is by invoking an 
enhancement in the carbon abundance by a factor ~ 1 0. Observations of circumstellar absorption lines from C I and C II by FUSE, 
reported recently by Rob erge et alJ (120051) . seem to indicate that carbon may indeed be over-abundant by this factor. Whether the 
gas is produced at this abundance or has evolved to it is, however, not clear at the moment. In the following sections we proceed 
to explore additional braking mechanisms that act independently of the relative abundance pattern of the gas. 

3.2. Collisions with Charged Dust 

Dust grains immersed in a plasm a acquire a non-zero char ge ( Sn>itzerll9^J), a phenomenon observed in the interstellar medium 
(ISM) and the solar system (e.g., Mendis & Rosenbera ll994l lHoranvll99 6). We therefore expect grains in the (3 Pic system to 
be charged as well. This might contribute to slow down the ion fluid by enhancing the ion-grain collision frequency relative to 
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FIG. 3. — Contours of the ratio between the ion slow-down timescale and the radiative acceleration timescale (equations Hil and Il4l in the /3Pic disk (p and z 
are disk cylindrical coordinates). As an example, we adopt CII as field particles and Fell as test particles with injection velocity v = c t h.f (f = 1). In regions of 
the disk where metallic gas has been observed, the ion slow-down time is orders of magnitude shorter than the radiative acceleration time, indicating that ions are 
dynamically coupled to a single fluid. Results are similar when we adopt other elements as field particles. 



the case where grains are neutral. The mass in dust grains is expected to be larger than that of ionized particles, satisfying the 
requirement of high inertia. We therefore focus on the requirement of high collision frequency, which depends on the potential 
to which grains are charged to. 

The relative importance of different grain charging mechanisms is determined by environment conditions (e.g., 
iMendis &R osenbergl ll994l) . In all cases, an equilibrium electrostatic potential is reached by balancing positive and negative 
currents. Given the spectral type of /3Pic, and the fact that the disk is optically thin, the dominant source of positive charge 
is the ejection of photoelectrons by stellar UV radiation, while negative charge is provided by the collection of thermal elec- 
trons. The contribution of thermal ions and secondary electrons is negligible. As detailed in Appendix 151 we solve for th e 
equ ilibrium potential of grains in the /3Pic disk given our model spectrum, temperature profile, and electron density (S 32.1I & 
12. 2\ . Due to the exponential fall-off of photon flux in the UV wavelen gths , the potential depends weakly (in fact, logarithmi- 
cally) on all variab les except for the wo rk function W of the grain (eq. IB9> . More detailed calculations yield that silicate dust 
(W = 8 eV. IWein gartner & Draine 2001) is charged to a positive potential e<\> ~ 0.2 eV, while carbonaceous grains (W = 4.4 eV, 
IWeineartner & D raine 2001 ) achieve a much higher potential e</> ~ 2.5 eV, both values being fairly insensitive to the environment 
factors. High value for the latter group results from the fact that stellar photons with energy > 4.4 eV are much more numerous 
than those above 8 eV. In any case, both potentials are much higher than ~ k^T ~ 10~ 3 eV, the typical potential achieved when 
only thermal collection currents are included. Figure [5] shows the potential of silicates and carbonaceous grains normalized by 
ksT, as a function of disk cylindrical coordinates. The contribution of the ejected photoelectrons to the ambient electron density 
is negligible, hence we did not couple grain charging with the ionization balance calculation. 

Each dust grain is order s of magnitude more massive than a single ion, therefore grains can be considered to be immobile in 
space. As shown in 33.11 ions feel radiation force as an ensemble. They move with a velocity v relative to the grains, with an 
internal velocity dispersion c± (equation II 0i. In a frame in which the center of mass of ions is at rest, dust grains drift with 
velocity —v. This relative motion is dissipated in a dynamical friction timescale (equation[8}, with dust grains and ions as test and 
field particles, respectively. Since the frictional force felt by ions is equal and opposite to that felt by dust grains, the timescale 
for slowing down ions in the stellar frame is given by 

, f A i- • \ m ion k B Tv 

fid = t s (dust, ion) = , (16) 

ftiust 87rn dust Z^ ust Zf on e 4 i/(01nA 
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FIG. 4. — Contours of effective radiation pressure coefficient /3 c ff (equation ll5l , calculated by averaging the radiation force over all ionized particles, as a 
function of disk cylindrical coordinates p and z. Solar composition is assumed. Carbon, with /3 = 0, is important for reducing the value of /3 e ff. Increasing the 
ionization as one moves to larger distances (both outward and upward) accounts for the lowering of e g. However, /3 e ff > 0.5 for all regions where metallic gas 
has been observed, indicating that an external agent is required to brake a solar-composition ionic fluid. 
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FIG. 5. — Equilibrium electrostatic potential of spherical grains e</> (see Appendix |B]f or details on the calculation) normalized by k^T, as a function of disk 
cylindrical coordinates. Left panel shows silicates (W = 8 eV) and right panel carbonaceous grains (W = 4.4 eV). Due to strong photoelectric currents, both types 
of grains are charged to positive potentials much greater than ~ k#T, the value expected if only thermal collection currents are responsible for charging grains. 
Results shown here are robust and depend only logarithmically on electron density, stellar luminosity, and orbital separation. 



where gdust is the volume density of dust grains, gj on is that of ions, and the limit m t ^> ;«f was used in equation (|8j. The grain 
charge enters equation d!6i through Zd ust . The dynamical friction timescale (equation |8j takes only into account the effect of 
distant encounters, i.e., impact parameters much larger than the grain size, under the assumption that the cumu lative effect o f 
these collisions outweights that of close encounters (including impactive ones) in determining the drag f orce (Spitzer 1956). 
Also, this expression does not include collective effects among ions (e.g., Northrop & Birmingham 1990). Both these effects 
are negligible for the following reasons. First, impact parameters are of order b ~ Zd us t<? 2 /(ksT), thus most collisions are within 
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FIG. 6. — Function K-H(^) (equation ll7L which determines the sign of the net force acting on an ion. Arrows indicate the stable and unstable nature of the 
equilibrium velocities v,j r if t and v max , respectively. For any given value of K < //max, Vmax is the maximum velocity for which braking can occur, whereas v^jft is 
the drift velocity that ions reach in a time fy (equation ll6l . 



the range a^b <C Ad, where a is the grain size and Ad is the Debye screening length (consequently, we use b^ B = a ~ 1 /im 
in equationlTTV Second, we have uj p a/c t h ~ 10~ 3 , where lo p = {Airn e e 2 /m^) 1 ' 2 is the electron plasma frequency, thus collective 
effects are not likely to contribute significantly to the drag force. 
Under these assumptions, the equation of motion of a representative ion is given by 



dv m 
mi 0n — = (p e ff- 



dt 



1)- 



= (/3eff-D 



GMm u 
— 

GMm K 



fid 

K-H(0 



K 



where 



K = 



((3 ef f-l)GMm ion k B T 
8irr 2 n dust Zj ust Zl a e 4 lnA' 



(17) 



(18) 



The function //(£) is always positive, with a single maximum H m3X w 0.2 at £ sa 1, and asymptotic limits H(£) ~ 2£/(37r'/ 2 ) for 
1, and H(Q ~ l/(2£ 2 ) for £ 3> 1. Therefore, if K > H m - dx , no braking is possible, since the net force acting on an ion is 
always positive. In this case, as v increases, the relative contribution of the drag force decreases as 1 / v 2 , leading to a runaway 
solution. If K < // max , there exist two equilibrium velocities for which the net force on the ion vanishes, as shown in Figure [6] 
We denote them as v max and Vdrift (v max > Vdrift)- Ions with initial velocities Vdnft < v < v max will feel a drag force until they move 
with the velocity vdrift relative to the dust, while those with initial v > v max will be accelerated outward without bound. These 
behaviors reflect the nature of charged dust braking. 

We then calculate the value of K for the /3Pic dust disk, taking mj on = WFe, Z{ on = 1, and Zd us t<? 2 /a = e<f>. The term n^a 2 
is understood to be summed over all grain sizes. As this term also appears i n the dust optical cross section, its value can be 
obtained from fits to scattered and thermal emission dArtvmowicz et al.lll989l) . In the present work we use an updated fitting 
formula for the dust profile, which assumes a constant grain albedo 0.5, and is based on HST observations of Heap et al. (2000) 
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FIG. 7. — Constant K normalized by the maximum value of function H (equation ! 171 for 100% silicate (left) and 100% carbonaceous grains (right). Braking by 
charged dust requires K/H miiK < 1. While silicate grains are incapable of braking the ions, carbonaceous grains attain enough charge to exert a significant drag 
force and are thus an effective braking agent. 



(P. Artymowicz, private communication): 



/ 2\ r 
7T(fl ) Adust = — 

w 



PO 



n-i/2 



exp 



0.7 



(19) 



with po = 120AU, W = 6.6( j o// 9 o) ' 75 AU, and ro = 2 x 10~ 3 . The dust cross section (or differential optical depth) peaks at 
^120 AU, with an outward drop-off and an inner clearing region. 

The resulting value of K as a function of disk cylindrical coordinates is shown in FigureQfor 100% silicate (left) and 100% 
carbonaceous grains (right), normalized by i/ max . It can be seen immediately that silicates do not attain enough charge to make 
their interaction with ions significant. On the other hand, carbonaceous grains attain a much greater Coulomb cross section. 
Braking is most effective in the mid-plane region where dust density is the highest, the efficiency decreasing with increasing 
altitude above the midplane. However, at very high altitudes (z > 80 AU and p ^100 AU, roughly at the location where Call 
emission is detected), somewhat anti-intuitively, dust braking becomes effective again. This is explained by Figure|5] which shows 
increasing grain potential as altitude rises, a fact that is in turn explained by the exponentially decreasing electron density with 
altitude. If the disk is made of 100% cabonaceous grains, ions will be slowed down to velocities compatible with observational 
constraints, except very near the star (p < 50 AU) where the dust density is too low. Figure [8] shows the corresponding values of 
Vdnft and v max for 100% carbonaceous grains. 

Are dust grains expected to have a significant fraction o f carbo n? The answers seem to depend on the environment in which 
grains are found. For ISM grains, Weinga rtner & Draind J2001I) deduced a combination of carbon and silicates, with carbon 
taking up ^35% by volume, largely what one expects out of solar composition. Chondritic meteorites have only a few percent of 
carbon by volume, while most interpla netary dus t particles collected in the Earth stratosphere are much more carbon rich, with 
> 50% of carbon in some cases ( K eller et alJl 994). Similar enrichment is found in comet Halley. It is thought that chondrites are 
relatively poor in volatiles because they have undergone so me processing and are le ss primitive. Three pieces of work on (3 Pic 
may hold some clues to this question: albedo calculations dArtvmowicz et all 1989) suggest that dust in the outer disk is likely 
comprised of bright silicates (or ices) darkened by only a small amount (< 1%) of carbon materials; 5 silicate emission features 
are detected in the inner (3 Pic disk (< 20 AU, Weinberger et al. 2003), though absent from the outer part; and a super-solar 
abundance for carbon in the gas phase is indicated by a UV absorption study (Ro berge et al.l2005T) . 

What are the constraints one can impose on the relative abundance of carbonaceous grains if they are solely responsible for 
braking the metallic gas? Braking requires K/H max < 1 where K/H max oc l/«dust- From Figure[7] we have K/H max < 0.1 near the 
midplane region (z <10 AU) outside p > 100 AU for the case of 100% carbonaceous grains. So for this region, a carbon fraction 
~ 10% is sufficient for braking. This fraction rises with altitude, inversely proportional to K/H mdx . We recall again that the dust 
density profile is still very uncertain, therefore these numbers should be taken as a very rough estimate. 

3.3. Collisions with Neutral Gas 

Co llisions between neutral and ionized sp ecies have previously been investigated as a braking mechanism in the j3 Pic system 
(e.g. JBeust et alJl989HLagrange et alll998j). The dipole mo ment that ions induce on neutrals enhances the collision cross section 
relative to the neutral-neutral value (e.g. JBeust et al][l989l) . Previous studies, using neutral-neutral collisions and f3 values for 
neutral species, have estimated t hat ~50 M m of n eutral gas are required to slow down elements affected by radiation force to 
within observational constraints (Brandekeret al. 2004). Here we perform a new estimate on the required mass to account for 
observations, using the fact that ions are dynamically coupled and are the only species feeling the radiation force. 



5 However, this result is based on a spatially unresolved SED. Revisiting the problem using resolved images i Hea p et all 



I may yield new insights. 
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FIG . 8 . — Drift velocity Vdrift (left) and maximum velocity for braking v max (right) for 100% carbonaceous grains, as a function of disk cylindrical coordinates. 
The value of v^m falls below observational limits (~1 kms~') in the majority of the disk. 



Assuming that each ion-neutral collision absorbs m n v\ momentum from ions, wh ere m n is t he mass of neutral particles and Vj 
is the velocity of ions, the drag force exerted on an ion by neutrals is given by (e.g., Beust et alJl989h 



(20) 



where, again, the subscripts i and n refer to ions (test) and neutral (field) particles, respectively, v is the relative velocity between 
the species, and k is given by 



' 4a n Z;e 2 



(21) 



where a n is the polarizability of neutral particles. Values of a n for several neutral atoms and molecules are shown in Table [3] 
The equation of motion for a representative ion is then 



dv rn GMnii 

m-r = (Peff- 1) — 5 g n nnv. 

at r z 



This has a single stable equilibrium solution when 



V = Vdrift : 



(/3 eff -l)GMm, 
TTKg n r 2 



(22) 



(23) 



Assuming that during the time fdra? it takes to slow down the ion, the distance it travels is much less than the distance to the star 
(vdi-iftfdrag r), the time dependent solution to 1221 is given by (Liseau 2003) 

(24) 



V(f ) PS Vdrift + (V0 - Adrift) exp (-? /t dag ) , 

where vo is the initial relative velocity between the incident ion and the neutral gas, and 

;;?i 



fdr 



ag ' 



■KKQ n 



(25) 



is the characteristic timescale for momentum damping due to ion-neutral collisions. 

Requiring a fixed drift velocity everywhere in the disk and integrating the required mass density (from eauationl23l over a disk 
with vertical scale height H, we obtain a minimum mass for braking 



AGMrm f H ((3 eff - 1) 



KVdrift Jo 



-pdpdz 



= 3 x 10 _z M ffi I — 
m H 



, 1/2 / \ 1/2 / \ /a 1 l -1 

m„ \ ' / a m \ ' I rn; \ / O.lkms 



Vdrift 



(26) 



where for numerical evaluation we have taken H = 100 AU, pmm = 10 AU, and p max = 1 000 AU. This result is much lower than 
previous estimates, for two reasons. First, since we are braking the ion ensemble, we have /3 e ft ~5, as opposed to individual 
species with (3 ^300 (Nal, Cal), requiring therefore less neutral material. Second, as only ions need to be braked, the collision 
frequency is enhanced relative to the neutral-neutral case due to induced dipole moment on the neutrals, the increase being of 
the same order as the reduction in the (3 value, further reducing the amount of braking material required to satisfy observational 
constraints. However, we note that since equation J26b assumes a constant drift velocity everywhere in the disk, the implied 
density profile g n oc 1 /r 2 is unrealistic. Thus the value in equation (126b should be taken as a lower limit only. 
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TABLE 3 

POLARIZABILITIES FOR NEUTRAL ATOMS AND 
MOLECULES 



Substance c* n a n / m„ 

(1(T 24 cm 3 ) (10~ 24 cm 3 molg" 1 ) 



Experimental 11 



H 2 0.787 0.390 

H 2 1.501 0.083 

He 0.208 0.052 

CO 1.953 0.070 

Calculated b 

H 0.693 0.686 

C 1.755 0.146 

N 1.046 0.075 

O 0.678 0.042 



Uohnsod 12*553) 

b Calculated using density functional B3LYP and 
basis set aug-cc- pVDZ (Johnson 2005) 



Enforcing solar abundances everywhere in the disk, we can calculate the drift velocity using the density profile of 32.21 which 
implies a total disk mass of 0.1 M . Figure [9] shows the value of Vdrift as a function of disk cylindrical coordinates, assuming 
F e II as test particles and H I as field particles. In most of the disk, this drift velocity is well below the observational constraints 
of Brandeker et al. (2004). However, in the zone corresponding to high-altitude Call emission (p ~ 120 AU, z ~ 80 AU) we have 
Vdrift ~ 26kms _I , a factor 100 above what is observed. This discrepancy has two possible explanations: either our hydrogen 
profile falls off too steeply with height (due to, e.g., a wrong temperature), and/or the solar abundance assump tion is not valid. 
The l ow temperatures in the outer regions of the /3Pic disk allow for the existence of molecules (e.g.. iKamp & van Zadelhoff 
1200 ll) . A possible candidate for a braking medium could thus be water vapor sputtered from icy dust grains, which certainly 
does not need to follow solar abundance relative to the metallic gas. However, any braking medium other than hydrogen would 
require a higher total mass for a given drift velocity, since this element has the highest a n /m n ratio and thus highest k, as shown 
in Table|3] H2O, for instance, would need to be ^3 times more massive to achieve the drift velocities in Figure|9] 

Finally, we address the inertia requirement. Since the drag force on the ions is equal and opposite to the force exerted by the 
ions on the neutral gas, the requirement that the gravitation al fo rce on the latter is much larger than the force from the ions means 
(Asff- 1)« < On, or AsffMions < M min . Using the results of 5j2"2] we obtain M ions w 8 x 10" 4 M e , or M min w 40M ions ~ 10/3 e ffM ions , 
therefore the requirement of high inertia is satisfied even by the lower limit of equation i26\ . As mentioned before, any braking 
medium which follows a physical density profile is likely to have a mass higher than M m [ n . 

4. DISCUSSION 
4.1. Sensitivity to Gas Temperature 

Our above results are obtained by assuming a thermal equilibrium between gas and dust. How reasonable is this and how 
s ensitive are our results to this as sumption? 

Kamp & van Zadelhoff (2001) made a thorough study of the gas temperature in the /3Pic disk, taking into account relevant 
heating and cooling mechanisms. They were able to set an upper limit ^300 K for a solar composition gas, the major uncertainty 
being the unknown relative drift velocity between dust grains and gas particles, which they argue to be the dominant heating 
mechanism. Our study indicates that the gas disk can be significantly depleted in hydrogen relative to solar abundance. One 
naively expects that the gas will be cooler than the above upper limit, as cooling is dominated by line emissions from carbon and 
oxygen. This needs to be confirmed. 

The observed gas disk has a finite thickness, h/r ~0.2 at 100 AU (Brandeke r et al.ll2.004l) . Assuming that this has a thermal 
origin (as opposed to a dynamical origin), and adopting a constant h/r ratio for the whole disk, we obtain a temperature profile 



r ~ 9xl ° 4 H^A^J {-t)' (27) 

where m is the mean mass of a gas particle, and n%c that of a carbon atom. This yields T = 900 K at r = 100 AU, much ho tter than 
the T = 60 K value in the r gas = rd ust approximation, and lies well above the upper limit of Kamp & van Zadelhoffl J2001h . This is 
puzzling and deserves further study. 

Even adopting such a seemingly extreme temperature profile, we find that many of our conclusions remain unchanged. To 
reproduce the observed Nal density profile, the electron density at 100 AU rises slightly from n e = 3.7 cm" 3 to n e = 5 cm" 3 . The 
ionization balance depends on temperature weakly through the radiative recombination coefficients. 6 As a result, carbon, being 
intermediately ionized, is the only element that sees its ionization fraction changed significantly (rises from ^50% to 80-90%). 



m \ / h/r\ ( AU 



6 Given the low gas densities, collisional ionization, charge exchange recombinations and dielectronic recombination are irrelevant for the ionization balance. 
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FIG. 9. — Magnitude of equilibrium drift velocity obtained by using the hydrogen density profile of i|2.2l (for solar composition) and equation 1231 . This 
corresponds to a total gas mass of 0. 1 Mq . While metallic gas in the midplane region will drift slowly, satisfying observational constraints, the high-altitude 
regions where Ca II has been detected (p ~ 120 AU and z ~ 80 AU) is flowing outwards too fast. 



This leads to a small decrease in the value of /? en ~ (more CII ions which do not see radiation pressure). The ratio £sAacc> relevant 
for ion coupling, increases by about an order of magnitude: the linear dependence on k%T in the numerator of equation (jSJi is 
softened by the increase in the ionized particles (mostly CII) and the slow increase of In A. Thus, coupling of ions remains 
efficient. The grain charge is mildly increased, as electron capture into grains is less efficient. However, the higher impact 
velocity increases the ratio K/H mdx by a factor between 3 to 5. Dust braking becomes less effective, although carbonaceous dust 
remains capable of braking gas in the mid-plane. Braking by neutral gas is little affected by changes in temperature, with a weak 
indirect dependence via /3 e ff. 

4.2. Gas Origin and Lifetime 

If either ion-grain or ion-neutral collisions are ultimately responsible for the slowing down of ions in the j3 Pic system, there 
is a finite drift velocity to which particles are slowed down. An estimate of the lifetime of the metallic gas can be obtained by 
calculating the time it takes for a single gas particle to travel the size of the disk R at a speed Vdnft- For charged dust, assuming a 
100% carbonaceous dust disk, we have 



1/2 /60K\ 1/2 ( 0.1 \ / R 



whereas for ion-neutral collisions we have 

, m „ M „ ,o« (as) " 2 (s.) " 2 (an) (-*-) ( ) ,29, 

\nif ) \a u J \m t J \300AU J \3xl0- 2 M @ ) y V ; 

One key question, important for constraints on gas evolution, is where the observed gas originates. Clues may be found in the 
observational evidence: 

(a) The metals observed in absorption follow approximately solar abundance (lLagrange et alll995l) . 

(b) There is a strong upper limit on the H2 column density in the disk from FUSE observations jLecavelier des Etangs et all 
120011) that corresponds to a disk mass ~0. 1 M®. 
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(c) The metallic gas spatial distribution is very similar to that of the dust dOlofsson et alJ200lllBrandeker et all2004l) . 

(d) The m etals observed in absorption and emission are close to rest relative to the star (Lagrange et alll998tlBrandeker et alJ 
l200l . 

As to the origin of the observed metallic gas, we see three different possibilities: 

1 . The gas could be primordial, that is, a remnant from the initial cloud out of which the star formed. Equation J29I then 
implies a gas mass of 30 M m to preven t the ions from leaving the system on a timescale shorter than the lifetime of the 
system, f age sa 10 7 yr (Zuck erman et al.ll2001l) . If the disk gas consists predominantly of hydrogen, then it would be in 
disagreement with (b). There remains the possibility that the gas disk is dominated by elements other than H, such as 
oxygen or oxygen-bearing molecules, which could have avoided detection. From an evolutionary point of view, however, 
it may be hard to justify why there should be such a large reserv oir of oxygen left, and little hy drogen. Furthermore, the 
recent dynamical modeling of dust behavior in a gaseous disk by Thebault & Augereau (2005) excludes a gas disk more 
massive than ~0.4M ffi . 

2. The gas could be produced by bodies evaporating as they f all into the star, and then blown/diffused out by the radiation 
pressure to large distances in the disk (Lag range et aljH998l) . Braking gas could not be produced in such a scenario, since 
the gas would then just accumulate close to the release radius near the star. The hypothesis of remnant primordial gas 
braking the metallic ions from the inner disk faces a fine-tuning problem: with too little braking gas the metals will leave 
too quickly, whereas with too much braking gas the metals would never have reached the outer disk. 

The productio n rate of gas required, M = M m n /t]\f R ~ 10~ 8 - 10~ 9 M m yr" 1 , is con sistent with the production inferred from the 
FEB scenario (Lagrange -Henri et al.ll9 88). According to Thebault et al. (2003), this high evaporation rate is inconsistent 
with the structure of the inner dust disk around (3 Pic. A lower evaporation rate could still be consistent with the observed 
metallic gas density profile, provided that the drift velocity is decreased by the same amount by increasing the braking gas 
density. Constraint (b) limits the braking gas increase to a factor < 3. 

If charged dust is primarily responsible for the gas braking, then (c) is a natural consequence of this. However, braking by 
dust is unlikely in this case. The main reason is that the gas is initially injected with a velocity 3> v max (Figure|SJ, even in 
the favorable case of carbonaceous grains. 

3. The gas could be cre ated in grain-grain collisions, in line with the production of small dust grains from minor bodies 
(Thebault et al. 2003). The correlation between the dust and gas spatial distribution (c) would be a natural consequence of 
the production mechanism. Braking by dust is possible but not required, since even if the metallic ions are rapidly depleted 
initially, the neutral low-/? gas is left behind until the appropriate abundance pattern is reached for the ions to self -brake. If 
this is the case, then the gas volatiles are predicted to be somewhat enhanced relative to the metals, and at the same time 
deficient in H, since grains are expected to contain very little hydrogen. 

If grains are able to brake ions, not only would the ions be tracing out the space distribution of smaller grains, they should 
also reveal the kinematics of dust. The near-Keplerian profile of metallic ions in the (3 Pic disk would therefore indicate 
that the grains, even at 300 AU, are still orbiting the star, raising the issue of how these grains get to these remote locations, 
since radiative blow out is not an option. 

In summary, our results indicate that the hypothesis that the observed metallic gas is primordial can be excluded. We cannot 
rule out the possibility that the metallic gas is produced by falling evaporating bodies and then diffused outwards, but this requires 
coincidental circumstances that need a priori justification. We conclude that the gas is most likely produced by grain evaporation 
(possibly due to grain-grain collisions) and predict a gas composition where metals are depleted relative to the volatiles, and 
where the volatiles have the same relative abundance as in the grains. 

4.3. Extension to other systems 

Given that several debris disks have been discovered around main sequence stars with different spectral types, we may ask 
whether the braking mechanisms we have studied so far are also important in those systems. We can do the simple exercise of 
calculating the radiation force coefficients, the ionization state, the value of /3 e ff, and attempt to draw some qualitative conclu- 
sions based on what we have learned in the previous sections. Table|4]shows radiation force coefficients and neutral fractions at 
r = 100 AU in the midplane, as a function of spectral type, for C, Fe, Si, Mg, S, Ni, Ca, Al, Na, and K. The spec tra are high resolu- 
tion PHOENIX atmosphere models (P. Hauschildt, private communication) similar to those described in 32. II flux calibrated and 
rotationally broadened according to the following stars: HR4796A (A0V, 3M Q , 55L Q , vsinz = 152kms _I ), HD 107146 (G2V, 
1M Q , 1 L Q , 200km s" 1 ), TWHya (K8V, 0.6 M , 0.1 L , Skins" 1 ), and AUMic (M1V, 0.5 M Q , 0.05 L Q , Tkms" 1 ). The ioniza- 
tion state for each case was obtained under the same constraints as the /3Pic data (A5V, 1.75M Q , 11 L Q , 130 kms -1 ), which is 
shown for comparison. 

As expected, while (3 values decrease with later spectral type, neutral fractions increase. Spectral types G2, K8 and Ml have 
radiation force coefficients much smaller than 1, except for Nal and KI. These species are expected to be depleted in stars 
around K8, since they are mostly neutral, while still having (3 > 0.5. Late spectral types also have emission lines arising from 
chromospheric activity, which we do not account for in these calculations. This might increase the (3 values for the M1V type, 
making Na and K depletion more likely in these systems. 

For the A0V case, all species are mostly in the ionized phase, but with higher (3 values than in (3 Pic. Using the densities of 
ionized species obtained by changing the spectral type to A0V, one obtains /3 e ff ~ 15 at r = 100 AU. This result should be taken 
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TABLE 4 

Radiation force, neutral fractions", and velocities before ionization for different 

spectral types 



iype 


l^usntity 


p 


Fe 


Si 


Ma 
IVlg 


c 


Ni 


Ca 


Al 


Na 




AOV 


A 


3 


60 


17 


200 


12 


48 


580 


99 


420 


200 




Ai 


0.1 


21 


25 


40 


0.01 


0.5 


180 


20 










ni/n u + 


0.01 


3E-6 


2E-7 


6E-6 


5E-4 


5E-6 


5E-7 


1E-7 


4E-5 


3E-5 




v ion [kms 4 ] 


2 


5E-3 


1E-4 


0.04 


0.02 


6E-3 


0.01 


2E-4 


0.06 


0.1 


A5V 


A 


0.03 


27 


6 


74 


0.6 


26 


330 


53 


360 


200 




Al 


2E-3 


5 


9 


9 


9E-5 


0.07 


53 


0.4 










ni/n u+ 


0.6 


9E-4 


1E-4 


6E-4 


0.06 


5E-4 


4E-6 


9E-7 


3E-4 


2E-4 




v ion [kms -1 ] 




0.5 


0.02 


1 


2 


0.3 


0.03 


5E-4 


3 


0.5 


G2V 


A 


1E-5 


0.06 


4E-3 


0.2 


2E-5 


0.1 


3 


0.3 


12 


26 




Al 


1E-7 


7E-3 


0.01 


0.06 


3E-8 


2E-6 


0.8 


6E-5 










ni/n a + 


0.6 


0.9 


0.06 


2 


0.2 


2 


6E-3 


1E-3 


0.08 


0.03 




v ion [kms 4 ] 














0.2 




17 


7 


K8V 


A 


6E-9 


4E-3 


4E-5 


2E-3 


3E-8 


7E-3 


0.2 


0.02 


3 


11 




Al 


2E-11 


2E-4 


4E-5 


2E-3 


2E- 


6E-8 


0.06 


3E-8 










"i/«ii+ 


0.7 


1 


0.07 


3 


0.2 


3 


0.5 


0.1 


14 


13 




v ion [kms -1 ] 


















9E+3 


2E+3 


M1V 


A 


7E-10 


5E-4 


6E-6 


4E-4 


8E-9 


7E-4 


0.04 


2E-3 


0.6 


2.1 




Ai 


4E-12 


2E-5 


7E-6 


2E-4 


6E-10 


2E-8 


0.01 


3E-9 










m/mi+ 


0.7 


1 


0.07 


3 


0.2 


3 


0.6 


0.1 


15 


20 






















1E+4 


2E+3 



"Neutral fractions w ere c omputed at r = 100AU in the disk midplane, under the same constraints as the 
values determined in i|2.2l for /3 Pic, changing only the spectrum. 

b Corresponding to model spectra (see text for details), flux calibrated and rotationally broadened according 
to particular cases (see text). No chromospheric emission was included in the model. 



only as a rough reference, since ionization fractions depend on the densities of different species, which in turn depend on disk 
mass and composition. 

Also shown are velocities before ionization Vi on calculated as in equation Q, which have a physical meaning only when (3 > 0.5 
(when < 0.5, the species feels a gravitational attraction weaker by a factor [1-/3], but is still bound to the star). Again, Na and 
K see m to b e the problem, being likely decoupled from the rest of the ions in the G2 disk, in the same way as Be is decoupled in 
/3Pic ( tj2.2\ . Ions are likely to be very strongly coupled for the AO case. 

5. CONCLUSIONS 

Motivated by the apparent contradiction between the strong radiation force acting on the gas in the (3 Pictoris disk and spatially 
resolved observations showing it to be consistent with Keplerian rotation, we have explored different braking mechanisms that 
are likely to operate in the disk. We have found that: 

a) All species affected strongly by radiation force are heavily ionized. The velocities that the short-lived neutral particles 
achieve as a result of radiative acceleration, before becoming ionized, are consistent with observational constraints. 

b) Ions are dynamically coupled due to their high Coulomb collision frequency. They feel a single effective radiation force 
coefficient /3 e ff. For solar composition, this coefficient is ^5, in which case ions cannot brake by themselves. If carbon is 
over-abundant by a factor > 10 on the other hand, the ion fluid may indeed be self-braking. 

c) If a significant fraction of the dust in the disk is carbonaceous, the ion ensemble can be slowed down to drift velocities 
below the local sound speed by photoelectrically charged dust grains. This fraction is as low as 10% near the midplane and 
rises with altitude above the midplane. 

d) Ions can also be slowed down by collisions with neutral gas (e.g., H2, H2O). The amount of neutral material required to 
satisfy observational constraint s on drift velocities is < 1 Mm, su bstantially less than previous estimates. The UV absorption 
upper limit on H% is ~ O.IM© (Lecavelier des Etangs et al. 2001). 

e) Reasonable assumptions on gas drift velocities imply gas lifetimes of order 10 4 - 10 5 yr, much lower than the age of the 
system and which require a replenishment mechanism, in analogy with the dust. 

The required gas production rate 10~ 13 M© yr" 1 is comparable to the FEB scenario prediction (e.g.. Lagra nge-Henri et alJl 988). 
although we argue that gas is more likely to be produced by local dust grain evaporation. A primordial origin of the gas is ruled 
out, since for having a lifetime of ~10 7 yr it would require a substantial quantity of neutral material, which would have already 
been detected. 

Further gas observations in this and other systems would be able to test several predictions arising from this work, such as the 
absence of P and Be in /3Pic, and the relative importance of ion-grain and ion-neutral collisions for slowing down ions. A more 
detailed study of the dust composition would shed light on the feasibility of the dust braking scenario. 
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APPENDIX 



UNCERTAINTIES IN THE RADIATION FORCE COEFFICIENTS 
The errors in the radiation force coefficients being quoted in Table^were calculated using the following expression: 



5 JL\ =( 5 lX iv 



So 



Pi 



SA» 



A i0 J \F ) 



limb./ 



(Al) 



where 6/3 is the error quoted in Tabled /?, and <5A;o are the radiation force coefficient and error in Einstein coefficient correspond- 
ing to the z'-th transition, respectively, (SF/F)]^^ is the fractional error in the stellar flux (at the given transition) due to uncertain 
limb darkening used in the rotational broadening of the spectrum, and (SF/F) ca \ is the fractional uncertainty in the overall flux 
calibration. The term (go/ ^2 go) represents the weighted average between different multiplets of ground states, when applicable 
(e.g., Fe I), where it is assumed that the population of a multiplet is proportional to its statistical weight. 

Values for (SA/A) range from 3% for the stronger transitions, to more than 50% for the weaker ones JMartin et al] ll999). 
The factor (SF/F)]^^ was estimated by rotating the atmosphere model with limb darkening coefficients e = {0,0.5,1} and 
calculating the fractional deviation, as a function of wavelength, relative to the e = 0.5 case. As mentioned in the main text, the 
error (SF /F) ca \ was estimated to be 4%, although systematic deviations may arise between the generic model spectrum and the 
particular ft Pic case. 

GRAIN CHARGING IN THE PIC DISK 
The relative importance of different grain charging mechanisms is determined by four parameters: the dust grain size a, 

the Debye screening length Ad of electrons, the mean intergrain separation re d , and the UV flux incident on the system 
(Mend is & Rosenberg! 19941) . A lower limit on the intergrain separation can be estimated by using the smalle st grain size and the 
highest dust mass in the disk. Assuming a higher limit for dust mass ~ 1 20 M ffi in the disk ( Artvmowicz l 19971) . using the smallest 
grain size a ~1 ^m (Backman & Paresce 1993), grain mass density ~1 gem" 3 , and a disk volume 7r(100AU) 2 x (10AU) we 

estimate « d '^ 3 > 180cm. Using the results of 32.21 we get Ad ~ 10-80 cm. Thus, the /3Pic disk is in the regime a < Ad < n^ 3 
where grai ns can be treated a s isolated particles immersed in an ionized gas, i.e., collective effects of charged grains play no role 
JMendis & Rosenberg! 19941) . 

The equilibrium electrostatic potential of a grain is found by balancing all the currents that charge the grain. Given the UV 
flux of the star, photoelectric charging is the dominant process, with a current per unit area given by (e.g, Draine 1978) 



Jph '■ 



(e<fi+W)/h 



hv-W 



f(E,hv)dE 



e(f> 



Q abs (v)Y(hv)-^dv 
hv 



(Bl) 



for the case of spherical grains with <ft > 0. Here, F v is the stellar flux per unit frequency, Q a b s the grain absorption efficiency, Y 
the photoelectric yield, W the grain work function, and 



f(E,hu) ■ 



hv-W V hv-W 



hv-W 



(B2) 



is the ph otoelectr on kine tic energy distribution function, which we approximate as a parabolic function (e.g., 
IWeingartner & Drainel20 01 ). The integral over electron energy is then 



hv-W 



f(E,hv)dE= 1-3 



+ 2 



(B3) 



Y(hv) = - I 1 



(B4) 



. hv-W J \hv-W 

The photoelectric yield is approximated as ( Drain eU978l) 

W_ 

2 V" ku, 

The fact that the grains in the disk have a minimum size of order ~1 /im (Back man & Parescell 19931) allows the use of work 

are W = 8 eV for silicates and W = 4.4 for carbonaceous grains 
rains also justifies taking the absorption efficiency for UV radiation as 
as Q d bs(v) = 1 (e.g.. lGreenberg|1971l) . 

The positiv e photoelectric current in equation JB 1 1 is balanced by a thermal electron collection current per unit area (e.g., 
lSDitzerll941l) : 



functions corresponding to bulk matter. The values adopted 
JWeingartner & Draind 2001). The minimum size of the grains 



J e = -es e rie 



2irm e 



1 + - 



kBT e 



(B5) 
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where we use a sticking coefficient s e = 0.5 JPrainell978l) . Given the low temp erature of the gas, the inciden t electrons have ener- 
gies <C 1 eV , not enough to make secondary electron emission important (e.g.. lMendis & Rosenberg! 19941). The ion c ounterpart 
to equation dB5b is found to be negligible, since these particles are strongly repelled when e<fi > (e.g., Sp itzeJl941l) . Provided 
that relative drift velocities between d ust g rains and electro ns are smaller than the el ectron thermal velocity (i.e., v < 40km s" 1 ), 
it is not necessary to correct equation dB5b for this effect (Northrop & Birmingham 1996), since positive charge is not provided 
by ions but by photoelectrons. We set T e equal to the gas temperature T, equation 0. 

The equilibrium grain potential is found by solving J p h~Je = 0, the result being shown in Figure [5] normalized by the local 
value of k^T . As mentioned in the main text, the contribution of the ejected photoelectrons to the ambient electron density is 
very small. We estimate an upper limit to the density of ejected photoelectrons as 



^e,dust ^ 



k*T 



akftT 



' 0.2 cm <n e 



(B6) 



for e(f)/(k B T) = 1000, T = 60 K, n e = 4cm" 3 , and a = 1 /im. In the following, we provide a rough estimate for the equilibrium 
potential, and show that it depends strongly on the grain work function (W), but logarithmically on all other parameters, making 
our results quite robust. 
In the Wien regime, photon flux can be approximated as 



v exp 



hv 
~~kvf~* 



(B7) 



0.7 eV for /3Pic, r is the distance to the star, and /o = F l/ {vo)exp(hvo/[kBT ir ])/i/Q, with F v (vq) rj 3 x 
Hz" 1 and hv^ =4.4eV. We integrate the photoelectric current (equation IB 1 i taking only account of the ex- 
ponential drop-off in the flux (and disregarding factors of order unity outside the exponential factor), 



where k^T^ - 
10~ 9 erg cm" 2 s 



P h 



h hv 



exp 



<j)+W 
k B T, 



AU\ 2 / W 2 /t B 7; 



h 4 



exp 



k B T+ 



(B8) 



hi/Kiecfi+W 

where we have used the fact that W > etfi. Even if this may not be a good approximation, the current depends exp onentially on 
the ratio (e<fi + W^/k^T* and other factors outside the exponential matter little. We also simplify the equation dB5l by assuming 
e<fi ~ k^T* ^S> k^T . We equate the above expression to the photoelectric current, and solve for the equilibrium potential: 

2 



ecf) = k B TJn 



1 



8s e «e 



foW 2 k B T* 
h 4 



T 



-W 



(B9) 



This shows that the charging potential depends on all environment variables logarithmically, and is only affected strongly by the 
work function. Substituting in our fiducial values (r = 100 AU, n e ~ 4 cm" 3 , s e = 0.5, T = 60 K) for the j3 Pic disk, we obtain 
e<j> ~ 2.8 eV for silicate dust (W = 8 eV) and ecf> ~ 5.5 eV for carbon dust (W = 4.4 eV). Both values are much greater than the 
values one expects if photoelectric charging is absent (e<j> ~ kgT ~ 10~ 3 eV). We numerically calculate these potentials (Fig. [5} 
using a realistic stellar spectrum and found them to be somewhat smaller than the values obtained here. 
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